Can correlations drive a band insulator metallic? 
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We analyze the effects of the on-site Coulomb repulsion f/ on a band insulator using dynamical 
mean field theory (DMFT). We find the surprising result that the gap is suppressed to zero at a 
critical Ud and remains zero within a metallic phase. At a larger Uc2 there is a second transition 
from the metal to a Mott insulator, in which the gap increases with increasing U. These results 
are qualitatively different from Hartree-Fock theory which gives a monotonically decreasing but 
non-zero insulating gap for all finite U. 



I 

o 
o 



in 



in 
o 



o 
o 



X 



There are many examples of interactions opening a 
gap in a metallic state: the Mott metal-insulator tran- 
sition, as well as spin or charge density wave gaps on 
nested Fermi surfaces. However, the possibility of inter- 
actions suppressing the gap in a band insulator all the 
way down to zero seems, at first sight, to be counter- 
intuitive. In this Letter, we show that this does occur 
in a simple model treated within a well-tested approx- 
imation scheme, namely dynamical mean field theory 
(DMFT) XM- 

Specifically, we start with a simple tight-binding band 
insulator with two bands, one filled and the other unfilled, 
and turn on an on-site Coulomb repulsion, the Hubbard 
U. We treat the effect of U using the DMFT within 
the paramagnetic sector, and show that with increasing 
U the gap of the band insulator first collapses to zero 
at a Uci and remains zero, leading to a metallic phase, 
for a finite range of U. This may be intuitively under- 
stood in terms of correlations screening out the one-body 
potential responsible for the gap in the band insulator. 
Eventually, with a further increase in U , there is a sec- 
ond phase transition at a Uc2 to a Mott insulator. Our 
calculated phase diagram is shown in Fig. ^ 

The model we consider has tight-binding electrons on 
a bipartite lattice (sublattices A and B) described by 



H 



(1) 



where t is the nearest neighbor hopping, U the Hub- 
bard repulsion and A a one-body potential which dou- 
bles the unit cell. The chemical potential is chosen 
to be ^ = ?7/2, so that the average occupancy is 
{{ua) + {ub)) /2 = 1 ("half-filling"). The Hamiltonian 
is sometimes called the "ionic Hubbard model" with 
A the "ionic" potential, and has been extensively studied 
especially in one dimension 0, 0, S IM S ^| j although 
the existence of an intermediate metallic state between 



the band and Mott insulators is apparently not settled 
conclusively even in the ID case. In any case, both our 
motivation and methodology are very different from pre- 
vious studies. 

We begin by considering various limiting cases at half 
filling. In the non-interacting [U = 0) limit the system is 
a band insulator with a band gap of A. However, a very 
different picture emerges in the atomic limit (t = 0). Here 
we obtain a "band insulator" for U /2 < A with site occu- 
pancies = and ub — "2 and a gap A — U/2 for either 
injecting or extracting an electron. This gap decreases 
with increasing U and collapses to zero at J7/2 = A when 
the singly-occupied states becomes degenerate with the 
empty one at site A and with the doubly-occupied one 
at site B. Beyond this point, for U/2 > A we get a 
"Mott insulator" with ha ^ hb = and a gap that 
goes like U when U ^ A. Hence, in this trivially soluble 
t = limit, we clearly see that the interaction U sup- 
presses the gap of (the state adiabatically connected to) 
the band insulator all the way down to zero, but only at a 
single point U = 2A. A question of considerable interest 
is whether this metallic point obtained in the atomic limit 
can broaden into a metallic phase when t is non-zero. 
Here we address this question using the DMFT ap- 

i is exact in the limit of 
been demonstrated to 
be very successful in understanding the metal-insulator 
transition P, Ql in the usual Hubbard model, which is 
the A = limit of eq. . We focus in this paper on the 
paramagnetic sector of eq. for which it is convenient 
to introduce the matrix Green's function 



proach. The DMFT approximation is 
large dimensionality and has be 



GQ^(k, iujn) 



CA(k,iuin) -Ck 
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where a, (3 are sublattice [A, B) indices, k belongs to 
the first Brillouin Zone (BZ) of one sublattice, icun — 
(2n + 1)ttT and T is the temperature. The kinetic 
energy is described by the dispersion and Ca{b) = 
iuJn T A + — S^(B)(*'^n)- Within the DMFT approach 
the self energy is purely local 1]. Thus the diagonal self- 
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energies Eq(iw„) are k-independent and the ofF-diagonal 
self-energies vanish (since the latter would couple the A 
and B sublattices). 

We note that within the Hartree-Fock (HF) approx- 
imation the self energy is = U {ua) /2 and the 
the difference in filling factors of the two sublattices 
Sn = {{nB) — {nA}) /2 is given by the T — Q self- 
consistent equation 5n = (A — U5n/2)^^^E^^ , where 
i?k = i/ej^ + (A - U5n/2Y is the HF excitation spec- 
trum. Thus, within HF theory there is a gap {A—U6n/2) 
for all A ^ and the system is a band insulator for all U. 
(For A = 0, one gets 6n = and the system is a metal 
for all U). However, a much richer, qualitatively differ- 
ent phase diagram emerges when we include fluctuations 
beyond HF, even at the level of the DMFT. 

The DMFT approach includes local quantum fluctua- 
tions by mapping P, Q the lattice problem onto a single- 
site or "impurity" with local interaction U hybridizing 
with a self-consistently determined bath as follows, (i) 
We start with a guess for Eq(w+) and 6n and compute 
the local GaiioJn) = Xlk ^""(^j *^") rewritten as 



GaiiuJn) = Ca(iw„) 



de- 



Po(e) 



(3) 



with a ~ A{B) and a = B{A) with po{e) is the bare 
density of states (DOS) for the lattice considered (see be- 
low). We actually need to solve the problem for only one 
sublattice and use the relations GAiitUn) = —GBi—itOn) 
and T,A{io-!n) = U — Sb(— ia;„) valid at half filling, (ii) 
We next determine the "host Green's function" [ij, |2| 
Qoa from the Dyson equation Goa{iuJn) = G~^{iu;n) + 
'E,a{ii-Un). (iii) We solve the impurity problem to obtain 
Tja{iujn) = Sq [Goaii^n)] (iv) We iterate steps (i), (ii) 
and (iii) till a self-consistent solution is obtained. 

We use as our "impurity solver" in step (iii) a general- 
ization of the iterated perturbation theory (IPT) 0, 0| 
scheme which has the merit of giving semi-analytical re- 
sults directly in the real frequency (ijj+ = uj + iO+) do- 
main. The IPT ansatz E^^^(w+) = + A„Ei^^(tj+) 
is constructed to be (a) exact for U/t <^ 1, (b) exact for 
t/U — 0, and (c) exact in the large u limit for all C//t, 
which imposes various exact sum rules. Here is the 
HF self energy with ria = —2 J°_^ Im Ga{u}^)duj/7^, and 
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(4) 

This has the form of the second order self-energy with 
Pa(ei) = -Im[goa(e+)]/7r, where ^o'J(t^+) = Goai^^) " 
E^-*^ is the Hartree corrected host Green's function and 

iV(ei,e2,e3) = /(ei)/(-e2)/(e3) + /(-ei)/(e2)/(-e3) 
where /(e) is the Fermi function. From condition (c) 
above we find that Aa = ^0(1 — na/2)/ [noa(l — noa/2)] 
with riQa — —2 j'^^lmQoaii^^)dLij/'K. Note that Aa is 




FIG. 1; The T — phase diagram of model of eq. (1) at half- 
filling obtained using DMFT for the Bethe lattice. The phase 
boundary Ud is denoted by open circles while Uc2 by filled 
circles. Inset: Detail of lower left hand corner of the phase 
diagram. 

same for both the sublattices. The results of this sim- 
ple approximate method are expected to be in semi- 
quantitative agreement with those of more exact but 
numerically intensive methods 0, 0, 0|- For simplic- 
ity, here we present the results for the T = solution 
of DMFT equations on a Bethe lattice of connectivity 
z — > 00. The hopping amplitude is rescaled as t ^ 
to get a non-trivial limit and the bare DOS is then given 
by Po(e) = \/4t^ — e^/ (27rt^) which greatly simphfies the 
integral in eq. (0) . We have also solved the DMFT equa- 
tions on the 2D square lattice and found that the results 
are qualitatively the samep^. 



1 0,4 



0.8 

1.4 


0.8 
1.4 




^ 0.8 Mott Insulator 
fo.4: 





U=7.5t 



-4 -2 2 4 6 
CO 



FIG. 2: The density of states p{u}) (in units of 1/t) plotted 
as a function of uj (in units of t). The three panels are all at 
A = O.lt for various values of U one for each phase. The gaps 
are clearly visible in the top (band insulator at [/ = O.lt) 
and bottom (Mott insulator at [/ = 7.5t) panels, while the 
middle panel for metallic state at U — 3.25t shows gapless 
excitations. 
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The phase diagram of Fig. ^ has been obtained from 
the Green's function, and the (renormahzed) density of 
states (DOS), given by p{uj) = - ImTr[G'(k, tj+)]/7r, 
calculated using the above procedure. Although the sub- 
lattice DOS Pa are not individually symmetric, = 
Pa{—u!) and thus p{oj) = + /c_b(w) is symmetric 

for half-filling. Fig. |3 shows how p{uj) evolves as a func- 
tion of U for a fixed A. At small U < Uci there is a 
gap in the spectrum, with the DOS at higher energies 
similar to the non-interacting (semicircular) result. We 
call this state a band insulator since it is adiabatically 
connected to the U = band insulator. At intermediate 
U the gap collapses to zero and we find a metallic 
state. The effects of correlations also show up in the pre- 
cursors to the upper and lower Hubbard bands at higher 
energies. At sufhciently large U > Uc2, the DOS again 
shows a gap and the system is a Mott insulator (adia- 
batically connected to the well studied Mott insulator at 
A = 0). The metallic phase, sandwiched between [/ci(A) 
and i7c2(A) in Fig.^ shrinks as A increases, consistent 
with there being a single metallic 'point' U — 2 A in the 
atomic limit. We return below to the question of the 
scale of the critical UcS as a function of A. 
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FIG. 3: Difference in filling factor 5n = (ub — nA)/2 for 
A = 0.025t plotted as a function of U/t: The DMFT result 
(filled circles) is much smaller than the HF result (triangles). 



The difference in filling 6n — (jib — nA)l1 and the en- 
ergy gap are plotted in Fig|31and Fig. 0] as functions of U 
for fixed A. For bn the DMFT results are much smaller 
than the HF results and become very small (of the or- 
der of error in evaluating (5n)for large values of U. But in 
the large U limit U ^ A, one can map the model (1) 
to the Heisenberg model with the spin exchange integral 
J = At^U/{U^ - 4A2) and show that 5n « IGt^A/C/^ 
is non zero. The DMFT result for the energy gap is 
qualitatively different from the simple HF result (mono- 
tonically decreasing, non-zero for all J7), and vanishes 
continuouslv|l3j at both metal-insulator transitions, at 
Uci and Uc2- 

Further insight into the scale of the gap and its non- 
monotonic variation with U is obtained by examining 
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FIG. 4: The energy gap for A = 0.025t plotted as a function 
of U/t. The DMFT gap (dots) from the DOS and Bgap (open 
squares) from eq. © are in excellent agreement. Note that 
the gap vanishes in the metallic phase at intermediate C/'s. 
(See text for discussion of error bars). In contrast, the HF 
results (triangles) do not show a metallic phase. 
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FIG. 5: Imaginary part of self energy T,'/^{u) plotted as a 
function of uj, where both are measured in units of t. All three 
panels are for A = Q.lt for various values of U. In the top 
(band insulator a.t U = Q.lt) and bottom (Mott insulator at 
U = 7.5t) panels E^(a;) vanishes identically for low energies 
\lo\ < 3i5gap. The metallic phase in the middle panel at ?7 = 
3.25t has a Fermi-liquid-like E^(a;) ~ w^. 



the self energy Sa(w+) = S^(tj) + ii;^(w). From Fig.H 
we see that S^(a;) vanishes for < 3i?gap in both 
the insulating phases. This can be understood from 
of the imaginary part of Sa of eq. I0J, which comes 
from a three fermion final state. In the metallic phase 
'S'^{uj) ^ Lo^ near = 0, characteristic of a Lan- 
dau Fermi liquid. In all the three phases, Ylj^iuS) can 
be written at low w as a Taylor expansion Sq(w) = 
S^(0) -I- (1 - Z-^) w + . . ., where Z can be shown to 
be independent of a. The spectral function is defined by 
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Aaa{^,^) = — l/7rImGac<(e, When S" is negligible 
we find, from eq. jSJ, that Aaa{£,oj) = 5(r(w) — e^) with 
r(cj) = + A-E^(u;))(w + ^ + A"E^(a;)). As e is 
real, w's which satisfy r{u!) < lie within the gap. The 
energy gap is then given by r (i?gap) = which, using 
the low-energy form of given above and the Kramers- 
Kronig relation, leads to the result 

-Bgap = Z|A- C/(5n/2 + 5| (5) 

where S — P dcuT,'^ {uj) / ttlo. We note that the energy 
gap obtained from the DOS (filled circles in Fig. ^ is 
in excellent agreement with that given by eq. (open 
squares in Fig. 0J, where the error bars on the latter 
are obtained from the estimated numerical error in E" at 
small UJ In particular we have checked very carefully 
that the gap indeed vanishes in the entire metallic phase 
within the limits of our numerical accuracy. 

We note that Eq. lO helps greatly in clarifying several 
features of our results(14l]. First, it shows exphcitly how 
correlations "screen" the one-body potential A responsi- 
ble for the gap in the band insulator and lead to a gap 
which is always less than the HF result in the band in- 
sulating phase. The dominant role in this being played 
by S{< 0) and by Z{< 1). Second, we estimate the 
location Ud of the band insulator to metal transition us- 
ing £;gap = which implies f/ = 2|A + S{U)\/5n{U). 
We find from our numerics that although |5'([/)| in- 
creases with [/, it is always much smaller than A. Thus 
Uci ^ 2A/^n([/ci) > 2A/5n(0), since 5n{U) decreases 
with increasing U (see Fig. 3). Thus for small A, we 
find that the very small 5n leads to Ud ^ A; e.g., for 
A = O.Olt we get Ud > lOOA ~ t, while for A = 0.025t, 
Ud > 20A ~ O.bt. Finally, the eventual increase in the 
gap with U in the Mott phase, can be also traced to the 
f7-dependences of Z and S in eq. 

We find that the low energy spectral function 
in the metallic phase A{e,uj) = — ImTr[G'(e, a;+)] ~ 
ZS{uj — Ze). Thus the Fermi Surface is given by e = 0, 
i.e., the same as for the metal with U — A = 0, a,s 
a consequence of the k-independence of DMFT self en- 
ergy. The quasiparticle residue Z, which is also the mass- 
renormalization, is a function of both U and A. Within 
the metallic phase it decreases with increasing U for a 
given A, and with increasing A for fixed U. 

In conclusion, we have analyzed a simple model of a 
band insulator with on-site Coulomb repulsion U. At the 
Hartree-Fock level the gap is reduced but the insulating 
behavior persists. However, when we treat correlation 
effects using DMFT we find the surprising result that the 
gap is suppressed to zero and there is a band insulator 
to metal transition at a critical value of Ud ■ At an even 
larger Uc2 there is a second transition from the metal to 
a Mott insulator. 

Our work raises several important questions. First, 
our results are obtained in the paramagnetic sector. We 



expect that the metallic phase will survive if antiferro- 
magnetism [l^ is suppressed due to frustration, in much 
the same way as in the DMFT treatment in the A = 
limit Q. It would be interesting to treat a model in 
which such frustration is explicitly included. Second, it 
is important to ask if in finite dimensions, i.e., with k- 
dependent self energies, other broken symmetry states 
might also appear, for example the bond ordered state 
which has been proposed as an intermediate phase in 
low dimensional studies 0, 0, 0| . Finally, it would be 
most interesting to look for experimental systems - either 
transition metal oxides [isj or even fermions in optical 
lattices where increasing correlations could drive a 
band-insulator metallic. All these are questions for fu- 
ture work. 

MR would like to thank N. Trivedi for very useful con- 
versations. 
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